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ABSTRACT 



Context. Low-energy cosmic rays are a fundamental source of ionization for molecular clouds, influencing their chemical, thermal 
and dynamical evolution. 

Aims. The purpose of this work is to explore the possibility that a low-energy component of cosmic-rays, not directly measurable 
from the Earth, can account for the discrepancy between the ionization rate measured in diffuse and dense interstellar clouds. 
Methods. We collect the most recent experimental and theoretical data on the cross sections for the production of and He^ by 
electron and proton impact, and we discuss the available constraints on the cosmic-ray fluxes in the local interstellar medium. Starting 
from different extrapolations at low energies of the demodulated cosmic-ray proton and electron spectra, we compute the propagated 
spectra in molecular clouds in the continuous slowing-down approximation taking into account all the relevant energy loss processes. 
Results. The theoretical value of the cosmic-ray ionization rate as a function of the column density of traversed matter is in agreement 
with the observational data only if either the flux of cosmic-ray electrons or of protons increases at low energies. The most successful 
models are characterized by a significant (or even dominant) contribution of the electron component to the ionization rate, in agreement 
with previous suggestions. However, the large spread of cosmic-ray ionization rates inferred from chemical models of molecular cloud 
cores remains to be explained. 

Conclusions. Available data combined with simple propagation models support the existence of a low-energy component (below ~ 
100 MeV) of cosmic-ray electrons or protons responsible for the ionization of molecular cloud cores and dense protostellar envelopes. 

Key words. ISM: cosmic rays, clouds - atomic and molecular processes 



1. Introduction 

Cosmic-rays (CRs) play a key role in the chemistry and dynam- 
ics of the interstellar medium (ISM). First, CR particles are a 
primary source of ionization, competing with stellar UV photons 
(absorbed in a thin layer of ~ 4 magnitudes of visual extinc- 
tion, McKee 1999) and X-rays produced by embedded young 
stellar objects (Krolik & Kallman 1983; Silk & Norman 1983). 
The ionization fraction in turn drives the chemistry of molecu- 
lar clouds and controls the coupling of the gas with the Galactic 
magnetic field (for a good review of the chemistry that occurs 
in the ISM in response to CR ionization see Dalgarno 2006). 
Second, CRs represent an important source of heating for molec- 
ular clouds because the energy of primary and secondary elec- 
trons produced by the ionization process is in large part con- 
verted into heat by inelastic collisions with ISM atoms and 
molecules. 

In general, the CR ionization rate in the interstellar gas de- 
pends on the relative amount of H, H2 and He (Dalgarno, Yan 
& Liu 1999). The first theoretical determination of the CR ion- 
ization rate was performed for clouds made only by atomic hy- 
drogen by Hayakawa, Nishimura & Takayanagi (1961). They 
assumed a proton specific intensity (hereafter, for simplicity, 
spectrum) proportional to the proton energy Ep for 0.1 MeV < 
Ep < 10 MeV and computed a; 4 x 10""' s"'. Spitzer & 
Tomasko (1968) determined a value (actually a lower limit) of 

> 6.8 X 10 s ' for HI clouds, assuming a CR proton 



spectrum declining below Ep ^ 50 MeV, and an upper limit 
of < 1.2 X 10"'^ s"', taking into account an additional flux 
of ~ 2 MeV protons produced by supernova explosions. To ob- 
tain the CR ionization rate of molecular hydrogen, a useful 
approximation is 1.5^^- ^ 2.3^^ (Glassgold & Langer 1974), 
giving ^ 10"'^ s"', in agreement with the lower limit on 
of Spitzer & Tomasko (1968). This value of is often refer- 
enced as the "standard" CR ionization rate in molecular clouds. 

A major problem in the determination of the CR ionization 
rate is that low-energy CRs are prevented from entering the he- 
liosphere by the solar wind and the interplanetary magnetic field 
(solar modulation). In practice. Earth-based measurements of 
CR fluxes give no information on the interstellar spectrum of 
protons and heavy nuclei for energies below ~ 1 GeV/nucIeon. 
Solar modulation also suppresses the flux of low-energy CR 
electrons, that shows considerable fluctuations akeady at ener- 
gies of 10-100 GeV (see e.g. Casadei & Bindi 2004). Since the 
cross section for ionization of molecular hydrogen by collisions 
with protons and electrons has a maximum at ~ 10 keV and 
~ 50 eV, respectively (see Sect. |2]l, it is clear that a knowledge 
of CR spectrum at low energies is an important limiting factor 
for an accurate calculation of the ionization rate in the ISM. A di- 
rect measurement of the shape of the CR spectrum at these ener- 
gies will be possible only when spacecrafts such as Pioneer and 
Voyager are well beyond the heliopause, the outermost boundary 
for solar modulation effects, believed to lie at 100-150 AU from 
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the Sun (at present, both Voyagers have already crossed the solar 
wind termination shock at 85-95 AU from the Sun). 

Over the last three decades, several values of rang- 
ing from a few 10"'^ s"' to a few 10"'* s"' have been 
obtained in diffuse interstellar clouds from measurements of 
the abundances of various chemical species, in particular OH 
(Black & Dalgarno 1977; Hartquist, Black, & Dalgarno 1978; 
Black, Hartquist, & Dalgarno 1978) and HD (van Dishoeck & 
Black 1986; Federman, Weber & Lambert 1996). However, the 
derived rates depend sensitively on several model assumptions, 
e.g. the value of specific chemical reaction rates and the intensity 
of the UV background. In dense molecular clouds, the determi- 
nation of the CR ionization rate is made even more uncertain 
by the sensitivity of molecular abundances to the level of deple- 
tion of the various species and the role of small and large grains 
in the chemical network. The values of derived by Caselli 
et al. (1998) in a sample of 23 molecular cloud cores (column 
density A?(H2) ~ lO^^ cm ^ ) through DCO+ and HCO+ abun- 
dance ratios span a range of about two orders of magnitudes from 
~ 10"'^ s"' to ~ 10"'^ s"', with a scatter that may in part reflect 
intrinsic variations of the CR flux from core to core. Finally, val- 
ues of of a few times 10"'^ s"' have been obtained in clouds 
of higher column density (A^(H2) ~ 10^^-10^'* cm"^) like the en- 
velopes surrounding massive protostellar sources (van der Tak & 
van Dishoeck 2000; Doty et al. 2002). 

The discovery of significant abundances of H^ in diffuse 
clouds (McCall et al. 1998), confirmed by follow-up detections 
(Geballe et al. 1999; McCall et al. 2003; Indriolo et al. 2007), 
has led to values of larger by about one order of magni- 
tude than both the "standard" rate and previous estimates based 
on the abundance of OH and HD in dense clouds. Given the 
relative simplicity of the chemistry of H^, it is now believed 
that diffuse clouds are characterized by CR ionization rates 
2 X 10""' s"' or larger This high value of i;^'- in the 
diffuse interstellar gas can be reconciled with the lower val- 
ues measured in cloud cores and massive protostellar envelopes 
by invoking various mechanisms of CR screening in molec- 
ular clouds due to either self-generated Alfven waves in the 
plasma (Skilling & Strong 1976; Hartquist, Doyle & Dalgarno; 
Padoan & Scalo 2005) or to magnetic mirror effects (Cesarsky & 
Volk 1978; Chandran 2000). An alternative explanation, based 
on the possible existence of a low-energy flux of CR particles, is 
that they can penetrate (and ionize) diffuse clouds but not dense 
clouds, as recently proposed by McCall et al. (2003; see also 
Takayanagi 1973 and Umebayashi & Nakano 1981). This latter 
scenario is explored quantitatively in the present paper. 

In this paper, we concentrate on molecular clouds, where 
hydrogen is present mostly in molecular form and we can ig- 
nore ionization of atomic hydrogen. In Sect. [8] we then apply 
our results to diffuse clouds, where the fraction of hydrogen in 
molecular form / = 2A^(H2)/[A^(H) + 2A^(H2)] has a mean value 
(/) ^ 0.6 (Indriolo et al. 2007), implying that the column den- 
sities of H and H2 are almost equal. This is justified because the 
quantity directly measured (or estimated) in the diffuse clouds 
examined in Sect. [8] is the ionization rate of H2 as derived from 
the measured abundance of Hj . 

The organization of the paper is the following. In Sect. |2l[3] 
and |4] we examine the ionization reactions of CR protons and 
electrons incident on H2 and He and other channels of electron 
production; in Sect. |5] we discuss the assumed interstellar spec- 
tra of CR protons and electrons; in Sect.|6]we discuss the energy 
loss mechanisms for CRs; in Sect. |7] we compute the ionization 
rate as function of the column density in a cloud; in Sect. [8] we 
compare our results with the available estimates of the CR ion- 



Table 1. CR reactions in molecular clouds 
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ization rate in diffuse and dense clouds; finally, in Sect. |9] we 
summarize our conclusions. 

2. CR reactions with I-I2 

CR particles (electrons, protons, and heavy nuclei) impact with 
atoms and molecules of the ISM producing ions and electrons. 
Table[T]lists the main CR ionization reactions involving H2 and 
He. In molecular clouds, a large majority of CR-H2 impacts 
leads to the formation of Hj via the ionization reaction 

^CR + H2 -> ^CR + H+ -i- e, (1) 

where kcR is a cosmic -ray particle of species k and energy Ei^, 
with cross section crj°" . Here we consider CR electrons (k - e), 
protons {k - p), and heavy nuclei of charge Ze (k = Z, with 
Z > 2). Low-energy CR protons, in addition, may react with 
ambient H2 by electron capture reactions, 

PCR + H2 -> H + H+, (2) 

with cross section cr^; . For an isotropic distribution of CR par- 
ticles, the production rate of Hj (per H2 molecule) is then 

= 47r y jk{Ek)[l + 0i(£t)]cr^™ (£^) Eu 

+4^ r '^"^ jpiE)(Ty {Ep)Ep, (3) 
Jo 

where jkiEk) is the number of CR particles of species k per unit 
area, time, solid angle and per energy interval (hereafter, we will 
refer to jkiEk) simply as the spectrum of particle k), /(H2) = 
15.603 eV is the ionization potential of H2, and fimax = 10 GeV 
is the maximum energy considered. The quantity 4>k{Ek) is a cor- 
rection factor accounting for the ionization of H2 by secondary 
electrons. In fact, secondary electrons are sufficiently energetic 
to induce further ionizations of H2 molecules, and their relatively 
short range justifies a local treatment of their ionizing effects. 
The number of secondary ionization produced per primary ion- 
ization of H2 by a particle k is determined by 

MEk) ^ -J— P"" P(Ek, E'Ja^riE:) E:, (4) 

where P{Ek,E'J is the probability that a secondary electron 
of energy E'^ is ejected in a primary ionization by a particle 
of energy Ek- The spectrum of secondary electrons declines 
rapidly with E'^ from the maximum at - (Glassgold & 
Langer 1973b; Cecchi-Pestellini & Aiello 1992). The function 
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4>e(Ee) giving the number of secondary ionizations after a single 
ionization by an electron of energy E,, has been computed by 
Glassgold & Langer (1973b) for energies of the incident elec- 
tron up to 10 keV. Above a few 100 eV, (p^ increases logarith- 
mically with Eg. For secondary electrons produced by impact of 
particles k, we adopt the scaling (f>k{Ek) ~ (peiEe = meEklnik) 
valid in the Bethe-Born approximation. Calculations by Cravens 
& Dalgarno (1978) confirm this scaling for protons in the range 
1-100 MeV. 

In the following subsections we summarize the available data 
for the ionization cross sections for proton and electron impact 
and for the electron capture cross section. The ionization of H2 
by CR heavy-nuclei (Z > 2) is computed in the Bethe-Born ap- 
proximation as described in AppendixlAl 

2.1. ionization of H2 by proton impact: 
PCR + H2 ^ pcR + H+ -H e 

The avalaible experimental data for proton-impact ionization of 
H2 have been summarized by Rudd et al. (1985). The cross sec- 
tion has a maximum at Ep x 70 keV and is considerably uncer- 
tain below ~ 1 keV. The data were fitted by Rudd et al. (1985) 
with expressions appropriate to the high- and low-energy re- 
gions, 

cr-=(crr'+cr,')-\ (5) 
where 

(Ti - ^na^Cx^ , (Th - 4TOg[A ln(l + x) + (6) 

with X = nieEp/mpIiH), /(H) = 13.598 eV, A = 0.71, B = 1.63, 
C = 0.51, D = 1.24. This expression is compared with exper- 
imental data in Fig.[T] For comparison, we also show in Fig.[T] 
the Bethe (1933) cross section for primary ionization of atomic 
hydrogen multiplied by a factor of 2. As it is evident, the Bethe 
formula reproduces very well the behavior of the ionization cross 
section at energies above a few tens of MeV. 
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10^ 10* 10= 10' 
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Fig. 1. Cross sections for proton impact on H2: ionization cross 
section cr™- (Rudd et al. 1985) and electron capture crj; ^- (Rudd 
et al. 1983) and total cross section cr^' for production of Hj. 
For comparison, the dot-dashed line shows the Bethe ioniza- 
tion cross section multiplied by a factor of 2. The two lower 
curves show the cross sections for dissociative ionization and 
double ionization of H2, multiplied by a factor of 10 and 100, 
respectively, obtained from the corresponding expressions for 
electron impact at equal velocity. Experimental data for the ion- 
ization cross section: stars, Gilbody & Hasted (1957); triangles, 
Afrosimov et al. (1958); diamonds. Hooper et al. {1961); filled 
circles, deHeer, Schutten & Moustafa (1966); Experimental 
data for the electron capture cross section: crosses, Curran, 
Donahue & Kasner (1959); empty circles, deHeer, Schutten & 
Moustafa (1966); asterisks, McClure (1966); squares, Toburen 
& Wilson (1972). 



2.2. Ionization of H2 by electron impact: 

ecR + H2 ^ ecR + + e 

The experimental data for electron-impact ionization of H2 
have been reviewed by Liu & Shemansky (2004). The absolute 
cross sections for electron-impact ionization of H2 measured by 
Straub et al. (1996) in the energy range E^ = 17 eV to - 
1 keV represent the currently recommended experimental val- 
ues (Lindsay & Mangan 2003). Analytic expressions and fitting 
formulae for the ionization cross section have been derived by 
Rudd (1991), Kim & Rudd (1994) and Liu & Shemansky (2004). 
Here we adopt the semi-empirical model by Rudd (1991) that 
gives an analytical expression valid up to relativistic velocities 
based on the theoretical results of Mott (1930), 



= 4-nalN 



/(H) 



/(H2) 



F(t)Git), 



where t - EelUVii), N - 2 (number of electrons of H2), 



n- 1 



1 +t 



nil 



1 - f l-«/2 

n-2 ' 



1 / M\ 
G(t) = -\^Ailnt + A2 + y-j. 



(7) 



(8) 



(9) 



with n = 2.4 + 0.2, Ai = 0.74 + 0.02, A2 = 0.87 + 0.05, 
A3 = -0.60 ± 0.05. For comparison, we also show in Fig. |2] 
the Bethe (1933) cross section for primary ionization of atomic 
hydrogen multiplied by a factor of 2. The Bethe formula repro- 
duces very well the behavior of the ionization cross section at 
energies above a few tens of keV. 



2.3. Electron capture ionization of H2: 
PCR + H2 ^ H + H+ 

In this charge-exchange process, a high-energy CR proton picks 
up an electron from the H2 molecule and emerges as a neutral H 
atom. The electron capture cross section has been fit by Rudd et 
al. (1983) with the expression 



Ana^AN 



/(H) 



/(H2) 



C + x'^+DxP' 



(10) 



where x — Ep/I(H), N-2 (number of electrons of H2), 
A = 1.044, B = 2.88, C = 0.016, D = 0.136, F = 5.86. This 
expression is compared in Fig. [T] with available experimental re- 
sults. 
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Table 2. Fit coefficients for the dissociative ionization and dou- 
ble ionization cross sections of H2 by electron impact. The cross 
sections are given by log(crp/10"'^ cm^) = 2^ an log(Ee/eV)". 



n 


a„ (diss, ion) 


a„ (doub. ion.) 





-53.23133 


-125.8689 


1 


96.57789 


172.0709 


2 


-67.57069 


-108.8777 


3 


23.32707 


34.18291 


4 


-4.004618 


-5.358045 
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0.272652 


0.335476 
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Fig. 2. Cross sections for electron impact on H2: ionization 
cross section o"™- (Rudd 1991), dissociative ionization o-diss.ion.^ 

and double ionization cross section cr'l°"^-^°"- (polynomial fits of 
Table|2] solidpait of the curves). For comparison, the dot-dashed 
line shows the Bethe ionization cross section multiplied by a fac- 
tor of 2. Experimental data for the ionization cross section: tri- 
angles, Rapp & Englander-Golden (1965); squares, Kossmann, 
Schwarzkopf & Schmidt (1990). Experimental data for the dis- 
sociative ionization cross section: diamonds, Straub et al. (1996). 
Experimental data for the double ionization cross section: filled 
circles, Kossmann, Schwarzkopf & Schmidt (1990). 



3. Additional reactions of CR electrons and protons 
with H2 

Additional ionization reactions that produce electrons are the 
dissociative ionization of H2, 



^CR + H2 ^ ^CR + H -H -H e. 



(11) 



with cross section crf^^- and the double ionization of H2, 



^CR + H2 ^ kcR + 2H+ + 2e. 



(12) 



with cross section cr^°*- These two processes contribute to 
the total CR production rate of electrons per H2 molecule. 



y f 



+4n y I jkiE)[l + MEk)]crf''- ™ (£*) d£, 

+8;r y r MEk)[l + Umcrf''''- '""(Et) d£,. (13) 

In the following subsection we examine the cross sections 
of these two processes for electron impact reactions, whereas 
for proton impact we assume cross sections equal to the corre- 
sponding cross sections for electrons of equal velocity. 



and 



^doub. io„.(£^) ^ ^doub. ion.(£^ ^ m,Ep/nip). 



(14) 



(15) 



As shown below, the cross sections of these processes are 
smaller by at least one order of magnitude than the correspond- 
ing ionization cross section, and the relative contribution of dis- 
sociative ionization and double ionization to the total electron 
production rate is expected to be small. 

3.1. Dissociative ionization of H2 by electron impact: 
ecR -I- H2 -> ecR -I- H -H -H e 



Absolute partial cross sections for dissociative ionization of H2 
by electron impact (threshold E'^"'^' '°" = 18.1 eV) have been 
measured by Straub et al. (1996) for incident electron energies 
ranging from E^ = 25 eV to E^ - I keV (see also Lindsay & 
Mangan 2003). Their results are in agreement with the reanal- 
ysis of Van Zyl & Stephen (1994) of the experimental results 
of Rapp, Englander-Golden & Briglia (1965), Krishnakumar & 
Srivastava (1994). For 1 keV < Eg <6 keV, the cross section has 
been measured by Takayanagi & Suzuki (1978). These measure- 
ments represent the currently recommended experimental values 
(Liu & Shemansky 2004). The data of Straub et al. (1996) and a 
polynomial fit of the data are shown in Fig. |2] The coefficients 
of the polynomial fit, valid for 18.1 eV < < 2 keV, are given 
in Table |2] 

3.2. Double ionization of H2 by electron impact: 

ecR + H2 — > ecR + 2H^ + 2e 

The energy threshold for this reaction is E'^°^'°- =51 eV. The 
cross section for this reaction is highly uncertain: the measure- 
ments by Edwards et al. (1988) and Kossmann, Schwarzkopf & 
Schmidt (1990) disagree by a factor of ~ 8. Here we adopt the 
latter set of measurements (shown in Fig.|2]i. The coefficients of 
a polynomial fit of these data, valid for 5 1 eV <£■(.< 4 ke V, are 
given in Table |2] 



4. CR reactions with He 

The CR production rate of He^ (per He atom) is 

^max 



f = An 



j,(£,)[l +<^i(£,)]cr'™(£,)d£, 



Jo 



jp{E)o--„'iEp)AEp 



(16) 



where /(He) = 24.587 eV is the ionization potential of He, crj°"- 
is the ionization cross sections of He for impact by particles k, 
and cr^ is the electron capture cross section. In the following 
subsections we describe the relevant cross section data proton 
and electron impact on He. As in the case of H2, the ionization 
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Fig. 3. Cross sections for proton impact on He: ionization cross 
section cr™- (Rudd et al. 1983), electron capture cross sec- 
tion crj; , and total cross section cr™ for production of He^. 
Experimental data for the ionization cross section: circles, Shah 
& Gilbody (1985). Data for the electron capture cross sec- 
tion: filled circles, Welsh et al. (1967); diamonds. Shah & 
Gilbody (1985). Data for the total ionization cross section: 
crosses, Pivovar & Levchenko (1967); triangles, Puckett & 
Martin (1970); squares, DuBois, Toburen & Rudd (1984). 



of He by CR heavy-nuclei is computed in the Bethe-Born ap- 
proximation described in AppendixlAl 

4. 1 . Ionization of He by proton impact: 
PCR -I- He ^ pcR + He+ + e 

Experimental measurements of He ionization by proton impact 
have been collected and fitted by Rudd et al. (1985). The cross 
section has a maximum at Ep x IQQ keV and is considerably un- 
certain below ~ 10 keV. Fig. [3] shows the available experimental 
data. We adopt the fitting formula of Rudd et al. (1985) given by 
eq. ^ and © with parameters A = 0.49, B = 0.62, C = 0.13, 
D = 1.52. 

4.2. Ionization of He by electron capture: 
PCR + He ^ H -h He+ 

The cross section for this charge transfer reaction has been mea- 
sured by Welsh et al. (1967) and Shah & Gilbody (1985). The 
cross section has a maximum at Ej, ^ 25 keV, where it is about 
one order of magnitude larger than the ionization cross section 
p-ion. ^ggg pjg Xotal ionization cross sections (cr™- + cr^- '^■) 
have been reported by DuBois, Toburen & Rudd (1984). 

4.3. Ionization of He by electron impact: 

ecR + He ecR + He^ + e 

Accurate experimental measurements of the cross section of ion- 
ization of He by electron impact are available (see Fig. |4]l and 
are in good agreement with theoretical calculations (Pindzola & 
Robicheaux 2000; Colgan et al. 2006). Here we adopt the fit- 
ting formula of Rudd (1991) given in eq. Q-® with N - 2, 
n = 2.4 + 0.3, Al = 0.85 + 0.04, A2 = 0.36 + 0.09, and 
A3 = -0.10 + 0.10. 




E [eV] 

Fig. 4. Cross section for He ionization by electron impact erf"- 
(Rudd 1991). Experimental data: triangles, Rapp & Englander- 
Golden (1965); circles Montague, HaiTison, & Smith (1984); di- 
amonds. Shah et al. (1988); squares, Kossmann, Schwarzkopf 8c 
Schmidt (1990). 



5. Local interstellar spectra 

From a theoretical point of view, if one assumes a uniform dis- 
tribution (in space and time) of CR sources characterized by a 
given "source spectrum" (usually a power-law in rigidity), CR 
propagation models can generate steady-state local interstellar 
(LIS) spectra resulting from a number of processes affecting the 
CR transport in the Galactic disk, like nuclear interactions, ion- 
ization energy loss, radioactive decay, escape from the Galaxy, 
etc. (see e.g. Berezinsky et al. 1990). These LIS spectra, in turn, 
can be used as input for solar modulation calculations to repro- 
duce the CR spectrum and the relative abundances of CR parti- 
cles measured at the Earth. The LIS spectra obtained in this way 
are clearly not uniquely defined, and a considerable range of LIS 
spectral shapes can be shown to be consistent with the measured 
CR flux with appropriate choices of parameters of the transport 
model (see e.g. Mewaldt et al. 2004, especially their Fig. 1). 

It is generally assumed that the LIS spectrum characterizes 
the energy distribution of CR everywhere in the Galactic disk, as 
long as the ISM properties do not depart from the uniform con- 
ditions assumed in the propagation model. With this assumption, 
Webber (1998) adopted LIS spectra for protons and heavy nuclei 
of energy greater than 10 MeV and electrons of energy greater 
than 2 MeV and combined them with data from Voyager and 
Pioneer spacecraft measurements out to 60 AU from the Sun 
to obtain a CR ionization rate x 3-4 x 10"'^ s"'. This is 
5-6 times the "standard" rate of Spitzer & Tomasko (1968) for 
atomic hydrogen. 

It is very uncertain, however, whether LIS spectra are really 
representative of the whole galactic disk, especially because the 
Solar System resides in a low-density (n ^ 10"^ cm"^) region 
produced by ~ 10 supernovae exploded over the past 10 Myr 
(the "Local Bubble"). In addition, to compute reliable CR ion- 
ization rates, the demodulated spectra need to be extrapolated 
down to ~ keV energies where the ionization cross sections have 
a maximum (see Sect.|2]|3]and|4li. Given these uncertainties, we 
discuss in the remainder of the paper the consequences for the 
CR ionization rate of making different assumptions about the 
low-energy behavior of CR spectra. In particular, we consider 
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Fig. 5. Proton LIS spectra of M02 and W98 (upper and lower 
solid curves, respectively). The dashed curves represent our 
power-law extrapolations of the spectra. For comparison, the 
cross sections for ionization of H2 by proton impact, electrons 
capture, and total ionization are also shown (in arbitrary units). 



for both protons and electrons a "minimum" and "maximum" 
LIS spectrum compatible with the available observational con- 
straints, and we compute the resulting ionization rates with the 
objective of comparing them with existing data for diffuse and 
dense clouds. 



5. 1 . Proton local Interstellar spectrum 

We consider two determinations of the proton LIS spectrum: 
Webber (1998, "minimum") and Moskalenko et al. (2002, "max- 
imum"), labeled respectively W98 and M02. Their characteris- 
tics are the following. 

(0 W98 estimated the LIS proton spectrum down to ~ 
10 MeV, starting from an injection spectrum parametrized as 
a power-law in rigidity, propagated according to the model of 
Webber (1987) and accounting for solar modulation following 
Potgieter (1995). The effects of solar modulation were refined 
using data from the Voyager and Pioneer spacecraft, then at dis- 
tances of ~ 60-70 AU from the Sun. The predicted LIS proton 
spectrum of W98 has a turnover around E x 100 MeV because 
of the dominant effect of ionization losses at low energies in the 
Galactic propagation model. Our extrapolation at low energies is 
a power-law in energy with exponent 0.95. 

(//) The adopted LIS spectrum of M02 (their "best-fitting" 
case) reproduces the observed spectrum of protons, antiprotons, 
alphas, the B/C ratio and the diffuse y-ray background. It is ob- 
tained from an injection spectrum which is a double power-law 
in rigidity with a steepening below 20 GeV, and a flattening of 
the diffusion coefficient below 4 GeV to match the B/C ratio 
at £ < 100 MeV. At low energies, our extrapolation follows a 
power-law in energy with exponent - 1 . 

Fig- H shows a comparison of the proton spectrum accord- 
ing to W98 and M02 (thick lines). The two spectra have been 
extrapolated as power laws down to ~ keV energies where the 
total ionization cross section, also shown in Fig. |5] has a broad 
maximum. 



5.2. Electron local Interstellar spectrum 

CR electrons (and positrons), although constituting a small per- 
centage of the corpuscular radiation, provide important informa- 
tion regarding interstellar propagation. This happens because CR 
electrons are more sensitive probes of ISM conditions than CR 
nuclei. In fact, electrons interact with: (/) the ISM, producing 
bremBtrahlung responsible for the largest part of galactic back- 
ground at y-frequencies; (//) radiation fields, generating radia- 
tion by inverse Compton scattering at X- and y-frequencies; (///) 
magnetic fields, producing synchrotron emission at radio fre- 
quencies. The electromagnetic radiation emitted by the interac- 
tion of CR electrons with other components of the ISM makes 
it possible to establish a relation between the observed radiation 
spectra and the energy distribution of the electrons. In particular, 
observations of the y-ray background in the 10 keV-100 MeV 
range, combined with measurements of the Galactic synchrotron 
spectral index in the frequency range 10 MHz- 10 GHz, provide 
indirect constraints on the LIS electron spectrum down to ener- 
gies of ~ 100 MeV. As for the proton spectrum, we extrapolate 
the LIS electron spectra to lower energies with power-laws to 
reach the peak of the ionization cross section at ~ 0.1 keV. Here 
we consider two different estimates of the LIS electron spectrum, 
both derived by Strong, Moskalenko & Reimer (2000). 

(/) The first spectrum, labeled COO, corresponds to the "con- 
ventional" model C of Strong et al. (2000), and is mostly derived 
from radio observations. It reproduces the spectrum of electrons, 
protons and alphas above ~ 10 GeV, satisfies the limits im- 
posed by positrons and antiprotons and the constraints on the 
synchrotron spectrum, but fails to account for the y-ray back- 
ground, especially for photon energies below ~ 30 MeV and 
above ~ 1 GeV. At low-energies, we have adopted a power-law 
dependence of the electron spectrum as Ef^^^. 

(ii) The second spectrum, labeled EOO, corresponds to model 
EOO of Strong et al. (2000). It reproduces the y observations at 
photon energies below ~ 30 MeV by a combination of bremB- 
trahlung and inverse Compton emission, assuming a steepening 
of the electron spectrum below ~ 200 MeV to compensate for 
the growth of ionization losses. The resulting increase in the syn- 
chrotron spectrum occurs at frequencies below 10 MHz, where 
the radio spectrum decreases abruptly due to the onset of free- 
free absorption. To fit OSSE data would require a LIS electron 
even steeper than EOO, but the excess y emission at ~ MeV en- 
ergies may be due to a population of unresolved point sources 
(Strong et al. 2000). At low energies, we have adopted a power- 
law extrapolation of the spectrum as E^^. 

In Fig. |6] we compare the two LIS electron spectra EOO and 
COO assumed in this work. 



5.3. CR Ionization rate for the local Interstellar spectra 

The values of ^f-, and ^^'^ per H2 molecule and He atom, 
respectively, obtained from numerical integration of eq. (O, ( fT3] l 
and ( fTSI l. with the jk{Ek) taken to be the adopted LIS spectra, are 
listed in Tab. [3] We have assumed a mixture of H2 and He with 
/h, - 0.83 and /hc = 0.17, corresponding to a He/H ratio of 0.1. 
We also list in Table |3]the energy density of each CR component, 
defined as 



Jo 



jk{Ek)Ek 
Vk(Ek) 



dEk 



(17) 



where jkiEk) is the particle's LIS spectrum and VkiEk) = 
c{E^/m^c'^ + 2Ek/nikC^y^^ /(I +Ek/mkC^) is the velocity of parti- 
cle k with kinetic energy Ek. We compute the total energy density 
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of CR as fi/t ~ (1 +^)&p, where ^ — 0.41 is the correction fac- 
tor for the abundance of He and heavy nuclei (see Appendix [All. 
The results listed in Table |3] suggest the following considera- 
tions: 

(/) Protons and heavy nuclei (plus secondary electrons) can 
produce ionization rates ranging from ~ 10"'^ s ' (in the case 
of the the spectrum W98, decreasing below Ep x 100 MeV) to 
~ 10"'"* S-' (spectrum M02, increasing below Ep ^ 100 MeV). 
The contribution of CR electrons to the ionization rate is negli- 
gible if the LIS electron spectrum flattens below E^ ~ 10 MeV 
(spectrum COO), but can become dominant if the spectrum in- 
creases at low energies. In practice, the ionization rate is propor- 
tional to the flux of CR particles in the energy range where the 
contribution to the integrals in eq. (O, ( fT3] ) and ( fTSI l is larger (see 
Sect.l7]andFig.[T4b. 

(//) The ratio of the CR ionization rate of He and H2 de- 
pends on the shape and absolute value of the assumed spectra. 
For CR protons, the ratio ^^''/^p ~ varies between 0. 15 (spectrum 
M02) and 0.64 (spectrum W98), whereas for electrons it varies 
between 0.38 (spectrum EOO) and 0.65 (spectrum COO). In gen- 
eral, since the ionization cross section for He decreases faster 
than that of H2 below the maximum, CR spectra rising with de- 
creasing energy result in a lower value of (^'^ j^^-. Given the 
sensitivity of modeled steady-state abundances of species like C, 
O2, H2O, H3 in dense clouds to the value of ^h<=/^H2 (Wakelam 
et al. 2006), it might be possible to constrain this ratio from a 
careful combination of molecular line observations and chemi- 
cal model predictions. 

(///) As anticipated, the CR production rate of electrons in 
molecular clouds is dominated by the CR ionization of H2 
(Sect.|2]l and He (Sect.|4]i. The contributions of dissociative ion- 
ization and double ionization to are small, about 5.5% and 
0.32% of the rate of production of electrons by single ionization 
of H2, respectively, independent of the adopted spectrum. 

(jv) The production rate of electrons, is generally larger 
than (but close to) the production rate of Hj . For the W98 pro- 
ton spectrum, the COO and EOO electron spectra, k 0.83- 
0.87^^. However, since we have included in the expression for 

the electron capture reaction (O whose cross section peaks 
at a lower energy than the ionization reaction ([T]i as shown in 
Fig. [U a CR proton spectrum rising at low energies may result 
in ^ > as in the case of the M02 spectrum. 

(v) With our assumed LIS spectra, the total CR energy den- 
sity varies from a minimum of 0.970 eV cm"-' (W98 plus COO) 
and a maximum of 1.80 eV cm"^* (M02 plus EOO), correspond- 
ing to an equipartition magnetic field of 6.2 //G and 8.5 jiG, re- 
spectively. These equipartition values are compatible with the 
"standard" value of the magnetic field of 6.0 + 1 .8 /iG in the cold 
neutral medium of the Galaxy (Heiles & Troland 2005). They 
have interesting consequences for the location of the solar wind 
termination shock (see discussion in Webber 1998). 

It is important to stress that the CR ionization rates listed 
in Table [3]have been obtained by integrating the spectra and the 
cross sections down to the ionization threshold of H2 and He, and 
they must therefore be considered as upper limits on the ioniza- 
tion rate. This is especially true for the electron spectrum EOO, 
which results in ionization rates exceeding the observed values 
by more than three orders of magnitude (see Sect. [8]). In the past, 
LIS spectra have been used to compute the CR ionization rate 
in the ISM assuming an appropriate lower cut-off in the CR en- 
ergy (e.g. Nath & Biermann 1994; Webber 1998). In this work, 
we use the LIS spectra to define the energy distribution of CR 
particles incident on the surface of the cloud. As we show in 



Table 3. CR ionization rates - and C^'^ (per H2 and per He, re- 
spectively), electron production rate f^, and energy densities &k 
of CR protons (p) and electrons (e) for the LIS spectra assumed 
in this work. The proton ionization rates include the contribution 
of heavy nuclei. 



k ref. <rr^ a 

(s"') (s"') (s"') (eVcm-^) 

"p W98 2.08 X 10-" 1.33 x 10"" 2.50 x 10"" 0953 

p M02 1.48 x 10"'^ 2.16x10-'-^ 3.49x10-'-^ 1.23 

~e COO 1.62 X 10-''' 1.05 x 10^''' \MxW" 00167 

e EOO 6.53x10-'^ 2.46 x lO^'^ 7.45 x lO^'^ 0.571 



10° 
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Fig. 6. Electron LIS spectra of EOO and COO (upper and lower 
solid curves, respectively). The dashed curves represent our ex- 
trapolations of the spectra. For comparison, the cross section for 
ionization section of H2 by electron impact is also shown (in 
arbitrary units). 



Sect.|6]and|2l the low-energy tail of the CR spectrum is strongly 
(and rapidly) modified by various energy loss processes when 
the particles propagate in a medium denser than the local ISM. 



6. Energy losses of CRs in the ISM 

The penetration of primary CR and secondary particles in in- 
terstellar clouds was studied by Takayanagi (1973) and more in 
detail by Umebayashi & Nakano (1981). In this paper we adopt 
the LIS spectra discussed in Sect. |5] to characterize the incident 
spectra and we follow the propagation of CR particles inside 
a molecular cloud with the so-called continuous-slowing-down 
approximation (hereafter CSDA) Q. In this approximation, the 
"degradation spectrum" of the CR component k resulting from 
the energy loss of the incident particles and the generation of 



' In the astrophysical literature this approximation is often referred 
to as the continuous energy loss regime, and when the propagation is 
dominated by these losses, it is often known as the thick target approx- 
imation (see e.g. Ramaty & Lingenfelter 1975; Ramaty, Kozlovsky & 
Lingenfelter 1996). 
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secondary particles is proportional to the inverse of the energy 
loss function, defined by 



Lk(Ek) 



1 



Im 

"n(H2) \ d£ 



(18) 



where n(H2) is the density of the medium in which the particles 
propagate and £ is the path length. Since we consider only energy 
losses in collisions with H2, our results are applicable only to 
clouds in which hydrogen is mostly in molecular form. 

In the following we consider CR propagation in molecular 
clouds assuming a plane-parallel geometry. It is convenient to 
introduce the column density of molecular hydrogen A^(H2), 

NiHi) = j n(H2)d{, (19) 

and to rewrite the energy loss function (eq. [TST i as 



Lk(Ek) 



dEk 



dA?(H2) 



(20) 



Let us then define jk(Ek,N) as the spectrum of CR particles of 
species k at depth A^(H2), with jk{Ek,Q) representing the LIS 
spectrum incident on the cloud's surface, defined by a column 
density A^(H2) = 0. To compute jk(Ek, N) we must consider all 
the processes that degrade the energy of the incident CR parti- 
cles. Assuming that the direction of propagation does not change 
significantly inside the cloud, it follows from eq. ( l20b that par- 
ticles of initial energy Ekx) reach energy Ek < Ekfi as a conse- 
quence of energy losses after propagating across a column den- 
sity A^(H2) given by 



N(n2) 



f 

Je, 



dEk 



■ n{H2)[Rk(Ek,Q)-Rk{Ek)l 



where Rk(Ek) is the range, defined as 

dEk 

at^ \ — 



rJ.) pE 

Rk(Ek)^ de^ 

Jei Jo 



-(dEk/dx) 



1 

^ Jo 



dEk 
LkiEkY 



(21) 



(22) 



Conservation of the the number of CR particles of each 
species implies 



jk(Ek,N)dEk^ j(Ekfl,0)dEk,o, 



(23) 



where, for a given value of A^(H2), the infinitesimal variation 
dEkfl of the particle's initial energy corresponds to an infinitesi- 
mal variation dEk of its energy at a depth A^(H2) given by 



dEk 



dE, 



k.O 



Lk(Ek) LkiEk.o) 



(24) 



(we ignore here that electron capture reactions of CR protons 
with H2 and He do not conserve the number of CR protons). 
Thus, the relation between the incident spectrum jk(Ekfl) and 
the spectrum at depth jkiEk, N) in the CSDA is 



jk{Ek,N)^jk{Ek,o) 



dEk 



dE, 



k.O 



- jkiEkfi) 



LkiEkfl) 
Lk(Ek) ■ 



(25) 



The energy loss functions for electrons and protons in H2 
are shown in Fig. |7] Some energy loss processes are com- 
mon to CR protons and electrons, like Coulomb interactions, 
inelastic collisions and ionization; others are peculiar to pro- 
tons (elastic collisions, pion production and spallation), others 
to electrons (bremfitrahlung, synchrotron emission and inverse 
Compton scattering). These processes are briefly reviewed in the 
following subsections. 



6. 1 . Energy loss of protons colliding with H2 

To determine the energy loss function of protons we have used 
the results collected by Phelps (1990) for energies in the range 
from 10"' eV to lO'' eV. For higher energies, between 1 keV and 
10 GeV, we have used data from the NIST Databas^l for atomic 
hydrogen multiplied by a factor of 2 to obtain the corresponding 
values for collisions with molecular hydrogen (NIST data do not 
include pion production at energies higher than about 0.5 GeV, 
that we computed following Schlickeiser 2002). The resulting 
energy loss function is shown in Fig.|7] The broad peak in Lp(Ep) 
at Ep X 10 eV is due to elastic collisions and to the excitation 
of rotational and vibrational levels, the peak at Ep « 100 keV to 
ionization, and the rapid increase at energies above ~ 1 GeV is 
due to pion production. For the low ionization levels character- 
istic of molecular clouds, the energy loss for Coulomb interac- 
tions of CRs with ambient electrons can be neglected at energies 
above ~ 1 eV (dashed line in Fig.|7]i. 

In Fig. [8] we show the quantity n(H2)RpiEp), obtained with 
a numerical integration of eq. ( l22l l. compared with data from 
the NIST Database at energies from 1 keV to 10 GeV. We also 
show the fit adopted by Takayanagi (1973) in a limited range of 
energies and the results of Cravens & Dalgarno (1978). As one 
can see, except for energies higher than ~ 100 MeV, where the 
NIST data do not include energy losses by pion production, the 
agreement between our results and the NIST data is very good. 



6.2. Energy loss of electrons colliding with H2 

To determine the electron energy loss function we have adopted 
the results of Dalgarno et al. (1999) for 10 - eV < < 1 keV 
and those of Cravens, Victor & Dalgarno (1975) for 1 eV < 
Eg < 10 keV. For higher energies, 10 keV < E^ < 10 GeV, we 
have adopted the loss function for electron-H collisions from the 
NIST Database multiplied by a factor of 2. The resulting energy 
loss function is also shown in Fig.|7] The first peak in Lf,(Eg) is 
due to the excitation of vibrational levels, the second to the exci- 
tation of the electronic levels and ionization, while at higher en- 
ergies the energy loss function is dominated by bremfitrahlung. 
As in the case of CR protons, we can neglect the contribution of 
Coulomb interactions for electrons at energies above ~ 1 eV. In 
Fig.[8l we show the range for electrons in H2, obtained as in the 
case of CR protons, compared with data from the NIST Database 
for lOkeV <E,<\ GeV. 



7. CR ionization rate in diffuse and dense clouds 

To compute the CR ionization rate in the ISM as a function of the 
column density A^(H2) of traversed matter, we follow the method 
of Takayanagi (1973). First, varying Ek and Ek.a from 0.1 eV to 
100 GeV, we determine the column density from the difference 
between RkiEk.o) and Rk{Ek). Second, tracing the level contours 
of the surface N{Ek.o, Ek) at different values of A^(H2), we ob- 
tain the relation between the energy of the incident CR particle, 
Ekfl, and the residual energy Ek, when the particle has covered a 
path inside the cloud corresponding to a given value of the col- 
umn density. We then fit the resulting Ek.o vs. Ek relation at fixed 
A^(H2) with the expression 



Ek.o(Ek,N) 



cEl + 



N 



Mb 



(26) 



|http : //physics . nist . gov/PhysRe£Data/Star/Text] 
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Fig. 7. Energy loss functions Lg(Ee) and Lp(Ep) for electrons 
and protons colliding with H2 (solid curves), compared with 
NIST data (circles); dashed curves show Coulomb losses for a 
fractional electron abundance n(e)/n(H2) - 10"^; dash-dotted 
curves labelled with n represent the energy loss by pion pro- 
duction computed following Schlickeiser (2002); dotted curves 
show the results by Phelps (1990) and Dalgarno et al. (1999) for 
P-H2 and e-H2, respectively. 




Fig. 8. Range Re(Ee) and Rp(Ep) for electrons and protons col- 
liding with H2 (solid curves), compared with NIST data (cir- 
cles) and the results of Cravens & Dalgarno (1978, squares); the 
dashed curve shows the fit by Takayanagi (1973) 



where Ek and Ekfi are in eV, and A^o in cm~^, b and c are non- 
dimensional. 

In Fig.|9][T0][TT] and [12] we show the CR spectrum obtained 
from eq. ( |25] l and ( |26] l for protons and electrons at values of 
A^(H2) ranging from 10'^ cm"- to 10^^ cm"-, inside a molecular 
cloud for the two incident spectra of protons and electrons de- 
scribed in Sect. |5] One can notice the correspondence between 
the shape of the proton spectra shown in Fig.|9]and[T0l and the 
energy loss function Lp(E) shown in Fig.|7] In fact, the relative 
minimum at about 10 eV in the attenuated spectrum corresponds 
to the energy loss peak due to elastic interactions and excitation 
of roto-vibrational levels, and the minimum at about 100 keV 



Fig. 9. Dashed curve: LIS proton spectrum W98 incident on 
the cloud's surface; dash-dotted curves: attenuated proton spec- 
tra at increasing depth in the cloud labeled by values of 
log[A?(H2)/cm-2]. 




10* 10° 

E [eV] 

Fig. 10. Dashed curve: LIS proton spectrum M02 incident 
on the cloud's surface; dash-dotted curves: attenuated proton 
spectra at increasing depth in the cloud labeled by values of 
log[A?(H2)/cm-2]. 



corresponds to the energy loss peak due to ionization. The same 
correspondence can be seen between electron spectra (Fig. [TTI 
and[T2]i and the energy loss function Le(Ee) (Fig. |7|l: the min- 
ima in the spectrum at about 1 eV and 100 eV are caused by 
the energy loss due to the excitation of vibrational levels, and 
to the excitation of electronic levels and ionization, respectively. 
This is a well-known property of the CSDA, where one approxi- 
mately obtains jk(Ek,N) oc l/Lk(Eit) independent on the column 
density if A^(H2) <sc A^o (see eq.l26]l. 

We are now able to calculate the CR ionization rate inside 
a molecular cloud as a function of the column density, with the 
attenuated spectra given by eq. ( |25] ). We compute the CR ioniza- 
tion rate for A^(H2) between 10'^ cm"^ and 10^^ cm"^, and we 
show the results for the four incident LIS spectra in Fig. [13] 

As a result of the detailed treatment of CR propagation, the 
decrease of the ionization rate with increasing penetration in the 
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Fig. 11. Dashed curve: LIS electron spectrum COO incident Fig. 13. CR ionization rate - as a function of the column den- 
on the cloud's surface; dash-dotted curves: attenuated proton sity A^(H2). Solid curves: contribution of CR protons (spectra 
spectra at increasing depth in the cloud labeled by values of W98 and M02); dashed curves, contribution of CR electrons 
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Fig. 12. Dashed curve: LIS electron spectrum EOO incident 
on the cloud's surface; dash-dotted curves: attenuated proton 
spectra at increasing depth in the cloud labeled by values of 
log[A?(H2)/cm-2]. 



Table 4. Fitting coefficients for eq. (|27T l. describing the atten- 
uation of the CR ionization rate for protons (p, also including 
heavy nuclei) and electrons (e). Eq. |27] | is valid in the column 
density range 10^° cm^^ < A^(H2) < 10"^ cm^^. 



k 


spectrum 


(s-') 




a 


P 


W98 


2.0 X 10- 


17 


0.021 


P 


M02 


6.8 X 10- 


16 


0.423 


e 


COO 


1.4 X 10- 


K 


0.040 


e 


EOO 


2.6 X 10- 


14 


0.805 



Table 5. Fitting coefficients for eq. (|28] | describing the attenua- 
tion of CR protons (p, also including heavy nuclei) and electrons 
(e). Eq.|28]l is valid for A^(H2) > lO^^ cm^^. 



k 


spectrum 












(s-') 




(g cm--) 


p 


W98 


3.4 X 10- 


li! 


44 


p 


M02 


5.4 X 10- 


18 


38 


e 


COO 


3.3 X 10- 


2(1 


71 


e 


EOO 


4.9 X 10- 


18 


35 



cloud at column densities in the range ~ 10^°- 10^^ cm"^ is char- 
acterized by a power-law behavior, rather than exponential at- 
tenuation, and can be approximated as 



A^(H2) 



1020 cm 



-2 



We have fitted this expression to the numerical results shown in 
Fig. [13] The coefficients ^o,* and a are given in Table HI The ex- 
ponential attenuation of the CR ionization rate sets in for column 
densities larger than ~ 10^^ cm-^, where depends essentially 
on the flux of CR particles in the high-energy tail of the incident 
spectrum (above ~ 0.1-1 GeV), and directly measurable on the 
Earth. In this regime, the attenuation of the CR ionization rate is 
expressed as function of the surface density of traversed matter 
2 = fimpN(li2), where nip is the proton mass and fi = 2.36 is the 
molecular weight for the assumed fractional abundances of H2 



and He (/h, =0.82 and /he 
fit the CR ionization rate as 



0.18). For S > 1 g cm ^, we can 



(27) /H^«fo,*exp- 



(28) 



where 2o ;(; is the attenuation surface density. In Table |5] we list 
the values of fo,/t and 'Lo.k obtained with the four spectra con- 
sidered in this work. The values for the attenuation surface den- 
sity listed in Table |5] are significantly lower than the "stan- 
dard" value of Nakano & Tademaru (1972) and Umebayashi & 
Nakano (1981), who obtain Eq ^ ^ 96 g cm"^ fQj- S > 50 g cm"^ 
(see also Umebayashi & Nakano 2009). 

It is important to stress that a large contribution to the ion- 
ization of H2 comes from low-energy protons and electrons con- 
stantly produced (in our steady-state model) by the slowing- 
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Table 6. CR ionization rate toward massive protostellar en- 
velopes 



Fig. 14. Differential contribution to the ionization rate £d^j ' / dE 
per logarithmic interval of kinetic energy, for the four spectra 
considered in this paper at a depth A^(H2) = 10^^ cm"^ [solid 
curves, protons; dashed curves, electrons). 



down of more energetic particles loosing energy by interaction 
with the ambient H2. In Fig. [14] we show the differential contri- 
bution of CRs protons and electrons to the ionization rate at a 
depth of A^(H2) - 10^^ cm"^, corresponding to the typical col- 
umn density of a dense cloud. For protons and heavy nuclei, the 
bulk of the ionization is provided by CR in the range 1 MeV- 
1 GeV and by a "shoulder" in the range 1-100 keV produced 
by slowed-down protons. This low-energy tail is produced dur- 
ing the propagation of CR protons in the cloud even when the 
incident spectrum is devoid of low-energy particles (as shown 
in Fig. |9] for the W98 spectrum). The largest contribution of 
CR-electrons to the ionization is distributed over energies in the 
range 10 keV-10 MeV, again reflecting the distribution of elec- 
trons in the propagated spectra (see Fig. [iTIandfTZb. Thus, the 
ionization rate at any depth in a cloud cannot be calculated by 
simply removing from the incident spectrum particles with ener- 
gies corresponding to ranges below the assumed depth. 

8. Comparison with observations 

To obtain the total CR ionization rate in molecular clouds, we 
sum the ionization rates of protons (corrected for heavy nu- 
clei as in Appendix lAb and electrons. With two possible spec- 
tra for each component, we obtain four possible profiles of . 
These are shown in Fig. [15] as function of A^(H2), compared 
with a compilation of empirical determinations of in dif- 
fuse and dense clouds. Our data sample includes: (/) diffuse 
clouds with A^(H2) from 10^*^ cm"^ to 10^^ cm"^ (14 detections 
and 15 upper limits, from Indriolo et al. 2007, including previ- 
ous data of McCall et al. 2002) and for the ( Per line-of-sight 
(Shaw et al. 2008); (//) molecular cloud cores with A^(H2) from 
10^' cm"^ to 10^^ cm"^ (data for low-mass cores from Caselli 
et al. 1998, WilHams et al. 1998, and for the prestellar core B68 
from Maret & Bergin 2007); (///) massive protostellar envelopes 
with A^(H2) from 10^^ cm"^ to 10^^ cm"^ (see Table|6]and refer- 
ences therein). 

The observational value of in diffuse clouds is obtained 
from the steady-state abundance of , produced by the CR ion- 
ization of H2 followed by a fast charge exchange reaction with 
H2 and destroyed mainly by electron recombination. The situ- 
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ation for dense molecular clouds and protostellar envelopes is 
more complicated. In the dense molecular gas, is removed 
by reactions with other molecules and atoms of the gas, e.g., by 
reaction with CO to form HCO^ and with O to form OH^. Thus 

can be determined from the measured abundance of a vari- 
ety of molecular ions such as HCO^, DCO^ and N2H^ (see e.g 
Caselli et al. 1998; van der Tak & van Dishoeck 2000; Doty et 
al. 2002). The resulting rates are nonetheless very uncertain, as 
they depend on the depletion of elemental C and O from their 
cosmic abundances, especially for clouds with a low degree of 
ionization, and are generally sensitive to the adopted chemical 
model. Here we adopt the values of derived by Caselli et 
al. (1998) with the data of Butner, Lada, & Loren (1995) and the 
chemical model of Leung, Herbst, & Huebner (1984). We view 
the range of values of obtained for different depletion factors 
as an indication of the associated uncertainties in the model de- 
terminations. In contrast with the study of Caselli et al. (1998), 
Williams et al. (1998) analyze molecular line data for a sample 
of low-mass cores using the chemical models of Bergin, Langer 
& Goldsmith (1995) and Bergin & Langer (1997). They con- 
clude that a single value (or a narrow range of values) of can 
reproduce reasonably well the observations for the majority of 
cores in their sample. 

The comparison between model results and observational 
data shown in Fig. [15] should be taken as indicative and inter- 
preted in a statistical sense, as also suggested by the large spread 
of values of at each value of A^(H2). First, the observational 
N(H2) is the entire column density through the cloud, whereas 
the model N{li2) is the column traversed by CRs incident over 
the cloud's surface. The exact relation between the quantities 
depend on factors like the cloud geometry and orientation with 
respect to the line-of-sight, and the variation of CR ionization 
rate with depth within the cloud. In addition, for the cloud cores 
of Caselli et al. (1998) we adopted the H2 column density esti- 
mated by Butner et al. (1995) from measurements of C'^O multi- 
plied by a factor of 2, to account for depletion of CO onto grains 
(Caselli et al. 1998). In fact, at the time of the study by CaselH 
et al. 1998, the almost complete disappearance of CO from the 
gas phase in cloud cores was still unknown. Second, many of 
the sight-lines where has been determined in diffuse clouds 
may have multiple cloud components, which would reduce the 
column density of a single cloud. It is probably safe to conclude 
that the observational column density is an upper limit to the 
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column density traversed by CRs incident on each cloud, and 
therefore the data shown in Fig. [15] should probably be shifted 
to the left by a factor or 2 or so. We will address the problems 
relative to cloud geometry and the effects of magnetic fields in 
a subsequent work. At any rate, from the comparison with ob- 
servational data, shown in Fig. [15] we can draw the following 
conclusions: 

(/) Although the gas column density of the object is by no 
means the only parameter controlling the CR ionization rate, the 
data suggest a decreasing trend of with increasing A^(H2), 
compatible with our models M02+C00, W98+E00, W98+C00. 
However, the measured values of are very uncertain, espe- 
cially in dense environments. Part of the large spread in the sam- 
ple of cloud cores may be due to a poor understanding of the 
chemistry. 

(//) The highest values of measured in diffuse clouds 
sight lines, could be explained if CR electrons are character- 
ized by a rising spectrum with decreasing energy. The BOO spec- 
trum represents an extreme example of this kind, and it results 
in values of somewhat in excess of the diffuse clouds ob- 
servations. The same spectrum accounts simultaneously for the 
CR ionization rates measured in most protostellar envelopes of 
much higher column density. Conversely, a spectrum of protons 
and heavy nuclei rising with decreasing energy, like the M02 
spectrum, can provide alone a reasonable lower limit for the CR 
ionization rate measured in diffuse clouds. 

(///) Without a significant low-energy (below ~ 100 MeV) 
component of electrons and/or protons and heavy nuclei, it is 
impossible to reproduce the large majority of observations. The 
combination of the COO spectrum for electrons with the W98 
spectrum for protons and heavy nuclei clearly fails over the en- 
tire range of column densities. Finally, a few molecular cloud 
cores and one dense envelope characterized by < 10"'^ s"' 
can only be explained by invoking the CR suppression mecha- 
nisms mentioned in Sect.[T]not considered in this work. 

In a recent paper, published after our work was completed, 
Indriolo, Fields & McCall (2009) analyze the implications of the 
CR ionization rate measured with in diffuse and dense clouds 
with an approach similar to that adopted in this paper. As in the 
present work, Indriolo et al. (2009) reach the conclusion that a 
low-energy CR component, likely produced by weak, localized 
shocks, can account for the relatively high CR ionization rate 
measured in diffuse clouds. However, their "best-fitting" CR- 
proton spectrum increases below Ep x IQO MeV as Ej^-^^, much 
more steeply than our steepest spectrum (M02, increasing as 
E^^). Indriolo et al. avoid a too large CR ionization (and heating) 
rate by cutting off the proton spectrum below Ep - 2 MeV and 
Ep - 10 MeV for diffuse and dense clouds, respectively. These 
energies correspond roughly to proton ranges in H2 of 10^' cm"^ 
and 10^^ cm"^, respectively. However, our work shows instead 
that when CR propagation is properly taken into account, low- 
energy particles are continuously produced by more energetic 
particles when they slow-down by interacting with the ambient 
medium (see Sect.O, making significantly to the ionization rate. 
In addition, we found that the contribution of CR-electrons to the 
ionization of H2, neglected by Indriolo et al. (2009), can be sig- 
nificant or even dominant over the contribution of CR-protons 
and heavy nuclei, without violating the available observational 
constraints. 

9. Conclusions 

The comparison between our models and the observational data 
available for diffuse clouds, dense cores and massive protostellar 
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Fig. 15. Total CR ionization rate a function of A^(H2) ac- 
cording to our models {solid curves). Observational data: filled 
circles and upper limits, diffuse clouds (Indriolo et al. 2007); 
square, diffuse cloud towards ^ Per (Shaw et al. 2008); empty 
circles, dense cores (Caselli et al. 1998); triangle, prestellar 
core B68 (Maret & Bergin 2007); filled squares, protostellar en- 
velopes (see Table |6]for references); (cross massive star-forming 
region DR21(OH) (Hezareh et al.2008). The box indicates the 
range of column densities and CR ionization rates compatible 
with the data analyzed by Williams et al. (1998). 



envelopes indicates that good agreement between theory and ob- 
servations can be obtained for the CR ionization rate of the ISM 
by including CR electrons with an energy spectrum increasing 
towards low energies, as also suggested by Webber (1998). In 
order to draw more stringent conclusions, it is necessary to use 
observational constraints derived from the ionization rates of dif- 
fuse clouds, as seen in Fig. [15] where the models differ mainly at 
low column densities. Our study points out the current limits to- 
wards a more accurate understanding of the ionization due to 
cosmic rays. There are essentially two limits: (1) the uncertainty 
in the CR spectrum at energies below ~ 1 GeV, and (2) the un- 
certainties in the empirically determined values of in diffuse 
and dense molecular clouds. The estimates of the ionization rate 
depend sensitively on the complex set of chemical reactions gov- 
erning the chemistry of the ISM, particularly on the so-called 
depletion processes that transfer molecules and ions from the 
gas phase to the solid phase. Despite these observational uncer- 
tainties, several important conclusions clearly emerge from our 
study: 

(0 values of measured in diffuse clouds are greater on 
average by an order of magnitude than those ones measured 
in dense molecular clouds. If confirmed, these data imply the 
presence of a CR proton and/or CR electron spectrum which in- 
creases at low energies. Thus, a combination of the spectra W98 
and COO for protons and electrons, respectively, is excluded by 
this set of observations; 

(//) values of measured in dense molecular clouds span a 
range of about two orders of magnitude and are subject to con- 
siderable uncertainty. It is difficult to establish how much of the 
observed spread is due to variations in the CR ionization rate. It 
is likely that in dense clouds the effects of magnetic fields on the 
propagation of CR particles cannot be neglected. In addition, it 
might be necessary to take into account the density distribution 
inside each cloud; 
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(///) the values of measured in massive protostellar en- 
velopes are somewhat higher than the predictions of our models 
at the corresponding column densities. This seems to suggest the 
presence of further ionization sources in these objects, as, for ex- 
ample. X-ray emission from the young stellar objects; 

(iv) The exponential attenuation of the CR ionization rate as- 
sumed in many studies is only established for column densities 
larger than ~ 10^^ cm"^. For the lower column densities consid- 
ered in this work, the ionization rate decreases as oc A^(H2)"" 
with a x 0.4-0.8 for the spectra that best reproduce the observa- 
tional data. 

It is reasonable to suppose that some of the observational 
uncertainties discussed here will be removed or reduced in the 
near future. With respect to the calculation of the ionization rate 
in dense clouds, the understanding of the complex chemical re- 
actions and of the depletion processes appear to be improving 
rapidly. Regarding the extrapolation of the measured CR spectra 
to low energies, the particle fluxes measured by Voyager and 
Pioneer spacecrafts outside of the solar magnetopause should 
provide important constraints on the energy distribution of CR 
protons and CR electrons, and, in any case, improve our under- 
standing of solar modulation. 
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Appendix A: Approximated corrections for 
lieavy-nuclei 



In the Bethe-Bom approximation, the cross section for the coUisional ionization 
of an atom or molecule depends only on the charge Zj. and the velocity of 
the incident particle. If A,t is the number of nucleons in the incident particle, the 
ionization cross section is (Tj(£j) a Z|(T(e), where e = EkIAk is the kinetic 
energy per nucleon and <T(e) is the same for all particle's species. If, in addition, 
the spectra of CR protons and heavy nuclei can be approximately desciibed by 
a single function such that jt(Et)AEt ~ /(iy(e)de, with representing 
the fractional abundance by number of species k, it is possible to reduce the 
calculation of the ionization rate by heavy-nuclei impact to that of protons, as 



jtiEt)[l+MEk)]<Tr {Ek)dEt ~ 



(!+'?) 



J;(H2) 



j„(Ep)[l+<l,p{E„)]cr'°''{E„)dE,„ 



where r/ is the correction factor for heavy nuclei ionization, 

7 *■ 

, jp 



(A.l) 



(A.2) 



Similarly, the correction factor to account for the energy density of heavy nuclei 
(eq. 117) is given by 



fp 



Ak. 



(A.3) 



Assuming for the CR abundance of heavy nuclei the standard solar abundance 
(Anders & Grevesse 1989), we obtain = 0.51, and f = 0.41, in agreement with 
the values = 0.50 and f = 0.42 of Indriolo et al. (2009). 
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